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Abstract. We have calculated long series expansions for self-avoiding walks and polygons on 
the honeycomb lattice, including series for metric properties such as mean-squared radius of 
gyration as well as series for moments of the area-distribution for polygons. Analysis of the 
series yields accurate estimates for the connective constant, critical exponents and amplitudes 
of honeycomb self-avoiding walks and polygons. The results from the numerical analysis agree 
to a high degree of accuracy with theoretical predictions for these quantities. 



1. Introduction 

Self-avoiding walks (SAWs) and polygons (SAPs) on regular lattices are among the most 
important and classic combinatorial problems in statistical mechanics. SAWs are often 
considered in the context of lattice models of polymers while SAPs are used to model vesicles. 
The fundamental problem is the calculation (up to translation) of the number of SAWs, c n , with 
n steps (SAPs, p n , of perimeter n). As for many interesting combinatorial problems, SAWs have 
exponential growth, c n ~ ^4^ n n 7_1 , where \i is the connective constant, 7 is a critical exponent, 
and A is a critical amplitude. A major challenge (short of an exact solution) is the calculation, 
or at least accurate estimation of, critical exponents and amplitudes. Here our focus is on 
the numerical estimation of such quantities from exact enumeration data. 

The success of series expansions as a numerical technique has relied crucially on several of 
Tony Guttmann's contributions to the field a asymptotic series analysis. In pioneering the 
method of differential approximants (see for a review and 'historical notes') Tony Guttmann 
has given us an invaluable tool which over the years has been proved to be by far the best 
(in terms both of accuracy and versatility) method for analysing series. In this paper we use 
long series expansions for self-avoiding polygons and walks on the honeycomb lattice to test 
the accuracy of various methods for series analysis. For the honeycomb lattice the connective 
constant, critical exponents and many amplitude ratios are known exactly, making it the perfect 
test-bed for series analysis techniques. 

The rest of the paper is organised as follows: In section |2] we give precise definitions of the 
models and the properties we investigate and summarise a number of exact results. Section 
contains a very brief introduction to the literature describing the algorithms used for the exact 
enumerations. In section |I] we give a brief introduction to the numerical technique of differential 
approximants and then proceed to analyse the SAP and SAW series clearly demonstrating 
how we can obtain very accurate estimates for the connective constant and critical exponents. 



Figure 1. Examples of a self-avoiding walk (left panel) and polygon (right panel) on the 
honeycomb lattice. 



Section [5] is concerned with the estimation of amplitudes. Not only do we obtain very accurate 
estimates for the amplitudes, but we also show how an analysis of the asymptotic behaviour of 
the series coefficients can be used to gain insight into corrections to scaling. Finally, in sectional 
we discuss and summarise our main results. 

2. Definitions and theoretical background 

An n-step self- avoiding walk u is a sequence of distinct vertices ... ,u n such that each 

vertex is a nearest neighbour of it predecessor. SAWs are considered distinct up to translations 
of the starting point ojq. We shall use the symbol f2 n to mean the set of all SAWs of length n. 
A self-avoiding polygon of length n is an n — 1-step SAW such that ujq and w n „i are nearest 
neighbours and a closed loop can be formed by inserting a single additional step between the 
two end-points of the walk. The two models are illustrated in figure ^ One is interested in the 
number of SAWs and SAPs of length n, various metric properties such as the radius of gyration, 
and for SAPs one can also ask about the area enclosed by the polygon. In this paper we study 
the following properties: 

(a) the number of n-step self-avoiding walks c n ; 

(b) the number of n-step self-avoiding polygons p n ; 

(c) the mean-square end-to-end distance of n-step SAWs (Rl) n ', 

(d) the mean-square radius of gyration of n-step SAWs (R 2 g ) n ', 

(e) the mean-square distance of a monomer from the end points of n-step SAWs (i?^) n ; 

(f) the mean-square radius of gyration of n-step SAPs (R 2 ) n ', and 

(g) the m th moment of the area of n-step SAPs (a m ) n . 

It is generally believed that the quantities listed above have the asymptotic forms as n — > oo: 
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The critical exponents are believed to be universal in that they only depend on the dimension of 
the underlying lattice. The connective constant fi and the critical amplitudes A-G^ vary from 
lattice to lattice. In two dimensions the critical exponents 7 = 43/32, a = 1/2 and v = 3/4 have 
been predicted exactly, though non-rigorously jSJIH]. In this work Nienhuis also predicted the 

exact value of the connective constant on the honeycomb lattice /i = + y/2. When analyzing 
the series data it is often convenient to use the associated generating functions such as 
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In the polygon generating function we take into account that SAPs have even length and the 
smallest one has perimeter 6. The SAW (SAP) generating function has a singularity at the 
critical point x = x c = 1/fi (x = x\ = l//x 2 ) with critical exponent —7 (2 — a). 
The metric properties for SAWs are defined by, 
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with a similar definition for the radius of gyration of SAPs. 

While the amplitudes are non-universal, there are many universal amplitude combinations. 
Any ratio of the metric SAW amplitudes, e.g. D/C and E/C, is expected to be universal [I]. 
Of particular interest is the linear combination |UE] (which we shall call the CSCPS relation) 
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where yt = 1/v and = 1 + 7/(2^). In two dimensions Cardy and Saleur [3] (as corrected by 
Caracciolo, Pelissetto and Sokal jSJ) have predicted, using conformal field theory, that H = 0. 
Cardy and Guttmann [0] proved that BF = aa®, where a is an integer constant such that 
p n is non-zero when n is divisible by a, so a = 2 for the honeycomb lattice, ao = 3y/3/4 is the 
area per lattice site on the honeycomb lattice. Richard, Guttmann and Jensen [Jj conjectured 
the exact form of the critical scaling function for self-avoiding polygons and consequently showed 
that the amplitude combinations G^ k 'B k ~ l are universal and predicted their exact values. The 
exact value for = i had previously been predicted by Cardy |Sj. 

The asymptotic form only explicitly gives the leading contribution. In general one would 
expect corrections to scaling so 
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In addition to "analytic" corrections to scaling of the form a,k/n k , where k is an integer, there 
are "non-analytic" corrections to scaling of the form 6fc/n Al+fc , where the correction-to-scaling 



exponent Ai isn't an integer. In fact one would expect a whole sequence of correction-to-scaling 
exponents Ai < A2 . . ., which are both universal and also independent of the observable, that is, 
the same for c n , p n , and so on. Furthermore, there should also be corrections with exponents such 
as nAj + mAj, etc., with n and m positive integers. At least two different theoretical predictions 
have been made for the exact value of the leading non-analytic correction-to-scaling exponent: 
Ai = 3/2 based on Coulomb-gas arguments 0111] and Aj = 11/16 based on conformal-invariance 
methods [§]. In a recent paper ^U] the amplitudes and the correction-to-scaling exponents for 
SAWs on the square and triangular lattices were studied in great detail. The analysis provided 
firm numerical evidence that Ai = 3/2 as predicted by Nienhuis. 

3. Enumerations 

The algorithm we used to enumerate SAPs on the honeycomb lattice is based on the finite-lattice 
method devised by Enting in his pioneering work, which contains a detailed description of the 
original approach for enumerating SAPs on the square lattice. A major enhancement, resulting 
in an exponentially more efficient algorithm, is described in some detail in ^2] while details of 
the changes required to enumerate area- moments and the radius of gyration can be found in • 
A very efficient parallel implementation is described in Jl]. The generalisation to enumerations 
of SAWs is straight forward as shown in j£] • An implementation of the basic SAP enumeration 
algorithm on the honeycomb lattice can be found in Most of the enhancements we made 
to the square lattice case can also be readily implemented on the honeycomb lattice. The only 
slightly tricky part is the calculation of metric properties (though the changes are very similar 
to those required for the triangular lattice |17j). 

Using the a parallel version of our honeycomb lattice algorithms we have counted the number 
of self-avoiding walks and polygons to length 105 and 158, respectively. For self-avoiding walks 
to length 96 we also calculate series for the metric properties of mean-square end-to-end distance, 
mean-square radius of gyration and the mean-square distance of a monomer from the end 
points. In fact the algorithm calculates the metric generating functions with coefficients c n (R 2 ) n , 
n 2 c n (Rg) n , and nc n (R 2 n ) n , respectively, the advantage being that these quantities are integer 
valued. For self-avoiding polygons to length 140 we calculate series for the mean-square radius 
of gyration and the first 10 moments of the area. Again we actually calculate the series with 
integer coefficients 8n 2 p n (R 2 ) n and p n (a k ) n - 

4. Differential approximants 

The majority of interesting models in statistical mechanics and combinatorics have generating 
functions with regular singular points such as those indicated in ((2j) and ©• The fundamental 
problem of series analysis is: Given a finite number of terms in the series expansion for a 
function F{x) what can one say about the singular behaviour which after all is a property of 
the infinite series. Without a doubt the best series analysis technique when it comes to locating 
singularities and estimating the associated critical exponents is differential approximants (see 
for a comprehensive review of differential approximants and other techniques for series analysis) . 
The basic idea is to approximate the function F{x) by solutions to differential equations with 
polynomial coefficients. The singular behaviour of such ODEs is much studied (see ^H]) and 
the singular points and exponents are easily calculated. 

A K th -ordev differential approximant (DA) to a function F(x) is formed by matching the 
coefficients in the polynomials Qi{x) and P(x) of degree iVj and L, respective, so that (one) of 
the formal solutions to the inhomogeneous differential equation 

^2Q i (x)(x±) i F(x)=P(x) 



Table 1. Critical point and exponent estimates for self-avoiding polygons. 
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0.29289321854(19) 
0.29289321875(21) 
0.29289321855(23) 
0.29289321859(19) 
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1.50000065(41) 
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1.50000038(33) 
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0.29289321861(37) 
0.29289321860(21) 


1.50000040(28) 
1.50000041(99) 
1.49999999(97) 
1.50000035(67) 
1.50000049(43) 



agrees with the first M = L + J2i(Ni + 1) series coefficients of F. Singularities of F(x) are 
approximated by the zeros xi of Qk(x) and the associated critical exponent A, is estimated 
from the indicial equation. If there is only a single root at this is just 

X = K-l - ( ^ K ~ 1 ^ 

The physical critical point is the first singularity on the positive real axis. 

In order to locate the singularities of the series in a systematic fashion we used the following 
procedure: We calculate all [L; Nq, N\, N2] and [L; Nq, Ni, N2, N3] second- and third-order 
inhomogeneous differential approximants with \Ni — Nj\ < 2, that is the degrees of the 
polynomials Qi differ by at most 2. In addition we demand that the total number of terms 
used by the DA is at least iV max — 10, where iV max is the total number of terms available in the 
series. Each approximant yields Nx possible singularities and associated exponents from the 
Nk zeroes of Qk( x ) (many of these are of course not actual singularities of the series but merely 
spurious zeros.) Next these zeroes are sorted into equivalence classes by the criterion that they 
lie at most a distance 2 k apart. An equivalence class is accepted as a singularity if it appears 
in more than 75% of the total number of approximants, and an estimate for the singularity 
and exponent is obtained by averaging over the included approximants (the spread among the 
approximants is also calculated). The calculation was then repeated for k — 1, k — 2, ... until a 
minimal value of 8 was reached. To avoid outputting well-converged singularities at every level, 
once an equivalence class has been accepted, the data used in the estimate is discarded, and the 
subsequent analysis is carried out on the remaining data only. One advantage of this method is 
that spurious outliers, some of which will almost always be present when so many approximants 
are generated, are discarded systematically and automatically. 

4-1- The polygon series 

First we apply our differential approximant analysis to the self-avoiding polygon generating 
function. In table ^ we have listed the estimates for the critical point x\ and exponent 
2 — a obtained from second- and third-order DAs. We note that all the estimates are in 
perfect agreement (surely a best case scenario) in that within 'error-bars' they take the same 
value. From this we arrive at the estimate x 2 c = 0.2928932186(5) and 2 - a = 1.5000004(10), 
where the error-bars reflect the spread among the estimates and the individual error-bars 
(note that DA estimates are not statistically independent so the final error-bars exceed the 
individual ones). The final estimates are in perfect agreement with the conjectured exact values 
x 2 c = l/fj? = 1/(2 + y/2) = 0.292893218813 ... and 2 - a = 3/2. 

Before proceeding we will consider possible sources of systematic errors. First and foremost 
the possibility that the estimates might display a systematic drift as the number of terms used 
is increased and secondly the possibility of numerical errors. The latter possibility is quickly 
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Figure 2. Plot of estimates from third order differential approximants for x 2 c vs. the highest 
order term used and 2 — a vs. x\. The straight lines are the exact predictions. 



dismissed. The calculations were performed using 128-bit real numbers. The estimates from a 
few approximants were compared to values obtained using MAPLE with 100 digits accuracy and 
this clearly showed that the program was numerically stable and rounding errors were negligible. 
In order to address the possibility of systematic drift and lack of convergence to the true critical 
values we refer to figure |5] (this is probably not really necessary in this case but we include the 
analysis here in order to present the general method) . In the left panel of figure |^1 we have 
plotted the estimates from third-order DAs for x 2 c vs. the highest order term iV used by the DA. 
Each dot in the figure is an estimate obtained from a specific approximant. As can be seen the 
estimates clearly settle down to the conjectured exact value (solid line) as N is increased and 
there is little to no evidence of any systematic drift at large N. One curious aspect though is 
the widening of the spread in the estimates around N = 140. We have no explanation for this 
behaviour but it could quite possibly be caused by just a few 'spurious' approximants. In the 
right panel we show the variation in the exponent estimates with the critical point estimates. We 
notice that the 'curve' traced out by the estimates pass through the intersection of the lines given 
by the exact values. We have not been able to determine the reason for the apparent branching 
into two parts. However, we note that the lower 'branch' contain many more approximants than 
the upper one. 

The differential approximant analysis can also be used to find possible non-physical 
singularities of the generating function. Averaging over the estimates from the DAs shows that 
there is an additional non-physical singularity on the negative x-axis at x = x_ = = 
—0.41230(2), where the associated critical exponent a_ has a value consistent with the exact 
value a_ = 3/2. In the left panel of figure 01 we have plotted a_ vs. the highest order term used 
by the DAs and we clearly see the convergence to a value consistent with a_ = 3/2. If we take 
this value as being exact we can get a refined estimate of x_ from the plot in the right panel of 
figure|Hl where we notice that the estimates for a_ cross the value 3/2 for x_ = —0.412305(5) 
which we take as our final estimate. From this we then get /i_ = 1.557366(10). 

4-2. The walk series 

Next we apply the differential approximant analysis to the self-avoiding walk generating function. 
In table |21 we have listed the estimates for the critical point x c and exponent 7 obtained from 
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Figure 3. Plot of estimates from third order differential approximants for the location x_ of 
the non-physical singularity and the associated exponent a_. The left shows a_ vs. the highest 
order term used and right panel a_ vs. X-. 



Table 2. Critical point and exponent estimates for self-avoiding walks. 
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1.34360(36) 
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1.3437685(54) 
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1.3437583(19) 
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1.3437669(53) 


0.54119610266(91) 


1.3437584(20) 


15 


0.5411961069(16) 


1.343776(68) 


0.5411961011(17) 


1.3437551(38) 
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0.5411961059(21) 


1.3437646(29) 


0.5411961022(26) 


1.3437580(59) 



second- and third-order DAs. Firstly, we note that estimates are about an order of magnitude 
less accurate than in the polygon case. Secondly, there are now small but nevertheless seemingly 
systematic differences between the second- and third order DAs (in particular the second-order 
homogeneous (L = 0) approximants are much less accurate than the other cases). On general 
theoretical grounds one would expect higher-order inhomogeneous approximants to be better in 
that they can accommodate more complicated functional behaviour. So based mainly on the 
third-order DAs we finally estimate that 0.541196102(4) and 7 = 1.343758(8). This is consistent 

with the exact values x c = l/fi= + y/2 = 0.541196100146 ... and 7 = 43/32 = 1.34375, 

though the central estimates for 7 are systematically a bit too high (and the second-order DAs 
are worse). 

In addition there is a singularity on the negative x-axis at x = —x c with a critical exponent 
consistent with the value 1/2, and a pair of singularities at x = ±0.64215(15)i with an exponent 
which is likely to equal 1/2 (note that the value 0.64215(15) is consistent with l//x_). These 
results help to at least partly explain why the walk series is more difficult to analyse than the 
polygon series. The walk series has more non-physical singularities and one of these (at x = —x c ) 
is closer to the origin than the non-physical singularity of the polygon series. Furthermore, 
as argued and confirmed numerically in the next section, the walk series has non-analytical 
corrections to scaling whereas the polygon series has only analytical corrections. All of these 
effects conspire to make the walk series much harder to analyse and it is indeed a great testament 
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Figure 4. Plot of estimates from third order differential approximants for 7 vs. the number 
of terms used by the DA and 7 vs. x c . The straight lines are the exact predictions. 



to the method of differential approximants that the analysis given above yields such accurate 
estimates despite all these complicating factors. 

In figure 0] we have plotted estimates for 7, obtained from third-order DAs, against the 
number of terms used by the DA (left panel) and against estimates for the critical point x c 
(right panel). The estimates for 7 display some rather curious and unexpected variations with 
the number of terms. Early on (around 80 terms) the estimates seems to settle at a value above 
the exact result. A little later the estimates start trending downwards so that around 95 terms 
they are in excellent agreement with the exact value. However, the estimates then unexpectedly 
start trending upwards again so that with more than a 100 terms the agreement with the exact 
result is only marginal. We also notice that in the right panel estimates for 7 vs. x c happen 
to just miss the intersection between the lines marking the exact predictions. This behaviour is 
curious and we have no ready explanation for it other than once again drawing attention to the 
quite complicated functional form of the generating function. The discrepancies between the 
estimates and exact values is marginal and certainly not significant enough to raise questions 
about the correctness of the theoretical predictions. 

4-3. Metric properties 

Finally, we briefly turn our attention to the series for the metric properties of SAPs and SAWs. 
We actually study the metric generating functions with integer coefficients 8n 2 p n {R 2 } n , c n (R 2 ) n , 
n 2 c n (R 2 ) n , and nc n {R 2 n ) n , which have critical exponents — (a — 1v) = —2, —(7+2^) = —91/32 = 
-2.84375, -(7 + 2zy + 2) = -155/32 = -4.84375, and -(7 + v + 1) = -123/32 = -3.84375, 
respectively (and as usual the polygon series use only the even terms). In table 01 we list the 
estimates obtained for the critical point and exponents using averages over third-order DAs. 
The exponent estimates from the SAP series are consistent with the expected value confirming 
v = 3/4 as are the estimates from the SAW series. The only possible exception is the end-to-end 
distance series where the estimates for both x c and the exponent are systematically a little to 
high. However, the discrepancy is not very large and probably not significant. 

5. Amplitude estimates 

Now that the exact values of \x and the exponents have been confirmed we turn our attention to 
the "fine structure" of the asymptotic form of the coefficients. In particular we are interested in 



Table 3. Critical point and exponent estimates for metric properties of SAPs and SAWs. 
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obtaining accurate estimates for the leading critical amplitudes such as A and B. Our method 
of analysis consists in fitting the coefficients to an assumed asymptotic form. Generally we must 
include a number of terms in order to account for the behaviour of the generating function at 
the physical singularity, the non-physical singularities as well as sub-dominant corrections to the 
leading order behaviour. As we hope to demonstrate, this method of analysis can not only yield 
accurate amplitude estimates, but it is often possible to clearly demonstrate which corrections 
to scaling are present. 

Before proceeding with the analysis we briefly consider the kind of terms which occur in the 
generating functions, and how they influence the asymptotic behaviour of the series coefficients. 
At the most basic level a function G(x) with a power-law singularity 

G(x) = ]T g n x n ~ A(x)(l - (6) 

n 

where A(x) is an analytic function at x = x c , gives rise to the asymptotic form of the coefficients 
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that is we get a dominant exponential growth given by fx, modified by a sub-dominant term 
given by the critical exponent followed by analytic corrections. The amplitude A is related to 
the function A{x) in © via the relation A = A(l/n) /T(£). If G(x) has a non-analytic correction 
to scaling such as 



G(x) = J2dnX n ~ (1 - fix)-* \A(x) + B{x){\ - /ix) A 
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we get the more complicated form 
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Figure 5. Plots of fits for the self- avoiding polygon amplitude B using in the left panel the 
asymptotic form (|12[) which ignores the singularity at x = X-, and in the middle panel the 
asymptotic form (|13j) which includes the singularity at x = x_. The right panel gives a closer 
look at the data form the middle panel. 
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Singularities in the complex plane are more complicated. However, a pair of singularities in 
the complex axis at ±i/r, that is a term of the form D(x)(l + t 2 x 2 )~ v , generally results in 

coefficients that change sign according to a + H pattern. This can be accommodated by 
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All of these possible contribution must then be put together in an assumed asymptotic 
expansion for the coefficients g n and we obtain estimates for the unknown amplitudes by directly 
fitting g n to the assumed form. That is we take a sub-sequence of terms {g n ,9n-i, ■ ■ ■ ,9n-k}, 
plug into the assumed form and solve the k + 1 linear equations to obtain estimates for the first 
few amplitudes. As we shall demonstrate below this allows us to probe the asymptotic form. 

5.1. Estimating the polygon amplitude B 

The asymptotic form of the coefficients p n of the generating function of square and triangular 
lattice SAPs has been studied in detail previously |19U12llTilll7| . There is now clear numerical 
evidence that the leading correction-to-scaling exponent for SAPs and SAWs is Ai = 3/2, as 
predicted by Nienhuis [2J 01 ■ As argued in ^3] this leading correction term combined with the 
2 — a = 3/2 term of the SAP generating function produces an analytic background term as can 
be seen from eq. (jHJ). Indeed in the previous analysis of SAPs there was no sign of non-analytic 
corrections-to-scaling to the generating function (a strong indirect argument that the leading 
correction-to-scaling exponent must be half- integer valued). At first we ignore the singularity 
at X- (since > x 2 it is exponentially suppressed) and obtain estimates for B by fitting p n 
to the form 



p n = fi n n-^ 2 



k 

B + J2 ai i ni 

i=l 



(12) 



That is we take a sub-sequence of terms {p n ,Pn-2, ■ ■ ■ ,Pn-2k} (n even), plug into the formula 
above and solve the k + 1 linear equations to obtain estimates for the amplitudes. It is then 



advantageous to plot estimates for the leading amplitude B against 1/n for several values of k. 
The results are plotted in the left panel of figure 03 Obviously the amplitude estimates are not 
well behaved and display clear parity effects. So clearly we can't just ignore the singularity at 
x_ (which gives rise to such effects) and we thus try fitting to the more general form 



Pn 



-5/2 



n 



+ (-i) n/ V 



n 



-5/2 



5>A 



(13) 



The results from these fits are shown in the middle panel of figure El Now we clearly have very 
well behaved estimates (note the significant change of scale along the y-axis from the left to 
the middle panel). In the right panel we take a more detailed look at the data and from the 
plot we estimate that B = 1.2719299(1). We notice that as more and more correction terms are 
added (k is increased) the plots of the amplitude estimates exhibits less curvature and the slope 
become less steep. This is very strong evidence that (|13[1 indeed is the correct asymptotic form 
of p n . 



5.2. Estimating the walk amplitude A 

From the differential approximant analysis we found that the walk generating function has 
non-physical singularities at x = —x c and x = ±i//i_. In addition we expect from Nienhuis's 
results (confirmed by extensive numerical work JHJ) a non-analytic correction-to-scaling term 
with exponent A = 3/2, and since 7 = 43/32 this correction term does not vanish in the walk 
case. Ignoring for the moment the pair of complex singularities we first try with to fit to the 
asymptotic form 



a.i/ri 



i/2 



+ (-1 



iv^~ 3/2 E 6 */ 



n 



(14) 



where the first sum starts at i = 2 because the leading correction is analytic. The resulting 
estimates for the leading amplitude A are plotted in the top left panel of figure El Clearly 
these amplitude estimates are not well behaved so we better not ignore the complex pair of 
singularities! Therefore we try again with the asymptotic form 



k k 

+ (-i) n //vr 3 / 2 ^&;/^ + (-i) [n/2l v-n- 3/2 Y, c i/ ni - ( 15 ) 

i=0 i=0 

The estimates for the leading amplitude A are plotted in the top right panel in figure El These 
amplitude estimates are not well behaved either and something is not quite right. Next we try to 
change the way we included the complex pair of singularities. In (|15[) we have assumed that all 
terms arising from the complex singularity have exactly the same sign-pattern. However, if we 
assume that the analytic correction terms arise from a functional form such as D(x)(l + /j?_x 2 ), 
where D{x) is a analytic, then the analytic correction terms would actually have a shifted sign- 
pattern. We therefore try fitting to the slightly modified asymptotic form 



k 
i=2 



i + Vo,/n 



i/2 



k k 

+ (-l) n /x n n- 3/2 bi/n* + ^ n~ 3 / 2 £(-1) L("+*)/ 2 J Ci /n\ (16) 
i=0 i=0 

where it should be noted that all we have done is change the way we include the terms from the 
complex singularities so as to shift the sign-pattern by a unit as i is increased. The new estimates 
for the leading amplitude are plotted in the bottom left panel of figure El and quite clearly the 
convergence is now very much improved. In the bottom right panel of figure El we show a much 
more detailed look at the data and from this plot we can estimate that A = 1.1449355(5). 
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Figure 6. Plots of fits for the self-avoiding walk amplitude A using different asymptotic forms. 
In the top left panel we show plots using the form l)14|) which ignores the complex singularity. 
In the top right panel we include the complex singularity via the form 1)15(1 while the bottom left 
panel includes the complex singularity via the alternative form 1)16(1. The bottom right panel is 
a more detailed look at this latter case. 



5.3. The correction-to-scaling exponent 

In this section we shall briefly show how the method of direct fitting can be used to differentiate 
between various possible values for the leading correction-to-scaling exponent Ai (recall the two 
theoretical prediction Ai = 3/2 by Nienhuis and Ai = 11/16 by Saleur). As already stated 
there is now firm evidence from previous work that the Nienhuis result is correct. Here we 
shall present further evidence. Different values for Ai leads to different assumed asymptotic 
forms for the coefficients. For the SAP series we argued that a value Ai =3/2 (or indeed any 
half-integer value) would result only in analytic corrections to the generating function and thus 
that p n asymptotically would be given by ()13|) . If on the other hand we have a generic value for 
Ai we would get 



\i n 



-5/2 



B + ^2a i /n i 



+ 



h/n 



Ai+i 



+ (-1 



,n/2 



fj,_n 



-5/2 



i=0 



(17) 



Fitting to this form we can then estimate the amplitude bo of the term 1/n 1 . We would expect 
that if we used a manifestly incorrect value for A\ then bo should vanish asymptotically thus 
demonstrating that this term is really absent from (|17|). So we tried fitting to this form using 
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Figure 7. Plots of estimates for the amplitude of the term 1/n 1 . The left panel shows 
results from fits to the form (|18|) where only the leading order term 1 /n Al is included (as well as 
analytical corrections). In the middle panel additional terms of the form l/n Al+ * are included 
and in the right panel terms like \/ n m ^ 1+l are included. 



the value Ax = 11/16. More precisely we fit to the generic form 

k k 
p n = ^ n - 5 / 2 £ ai/n a * + (-l) n/ V™ ™~ 5/2 E (18) 

j=0 i=0 

In the first instance we include only the leading term arising from Ai, that is we use the 
sequence of exponents a» = {0, 11/16, 1, 2, 3, . . .}. We also fit to a form in which the additional 
analytical corrections arising from Ai are included leading to the sequence of exponents 
ai = {0, 11/16, 1, 27/16, 2, 33/16, 3, 49/16, . . .}. As stated in Section |2] more generally one would 
also expect terms of the form 1 / n m Ai+« with m a non-negative integer. This leads to fits to the 
form given above but with a, = {0,0.6875,1,1.375,1.6875,2,2.0625,2.375,2.6875,2.75,3...}. 
The estimates of the amplitude of the term l/n Al as obtained from fits to these forms are 
shown in figure [7| As can be seen from the left panel, where we fit to the first case scenario, the 
amplitude clearly seems to converge to 0, which would indicate the absence of this term in the 
asymptotic expansion for p n . In the middle and right panels we show the results from fits to the 
more general forms. Again the estimates are consistent with the amplitude being 0. Though 
in this case the evidence is not quite as convincing. This is however not really surprising given 
that the incorrect value Ai = 11/16 gives rise to a plethora of absent terms which will tend to 
greatly obscure the true asymptotic behaviour. 

5.4- Amplitude ratios D/C and E/C 

From fits to the coefficients in the metric series we find AC = 1.0141(1), AD = 0.14225(5) and 
AE = 0.4458(1) and thus the ratios are 

D/C = 0.14027(6) and E/C = 0.43960(15) 

D/C and E/C can also be estimated directly from the relevant quotient sequence, e.g. 
r n = {Rg) n /{R 2 /) n , using the following method due to Owczarek et al. [2*U] : Given a sequence of 
the form g n ~ <7oo(l + °/ n + •■■)> we construct a new sequence {h n } defined by h n = Tlm=i 9m- 
The associated generating function then has the behaviour h n x n ~ (1 

— 9oo x ) 1 &nd we 

can now estimate form a differential approximant analysis. In this way, we obtained the 
estimates 

D/C = 0.1403001(2) and E/C = 0.439635(1) 
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Figure 8. Plots of the estimates for the amplitude combination 8BF. 

These amplitude estimates leads to a high precision confirmation of the CSCPS relation 
H = 0.000003(13). 

5.5. Amplitude combination BF 

Next we study the asymptotic form of the coefficients r n = 8n 2 p n (R 2 } n for the radius of gyration. 
The generating function has critical exponent — (a + 2u) = —2, so the leading correction-to- 
scaling term no longer becomes part of the analytic background term. We thus use the following 
asymptotic form: 



r n ~ M n 



8BF + ^2ai/n 



l+i/2 



i>0 



+ {-l) n ' 2 l i n _nY j b i /n l . 



(19) 



i=0 



In figure |S1 we plot the resulting estimates for the amplitude 8BF. The predicted exact value 
is BF = ^zo-a = = 0.0411312745..., where for the honeycomb lattice a = 2 and 

clq = 3v3/4. Clearly extrapolation of these numerical results yield estimates consistent with the 
theoretical prediction. 



5. 6. Amplitude ratios of area-weighted moments 

The amplitudes of the area- weighted moments were studied in • We fitted the coefficients to 
the assumed form 



np n (a m )n ~ ^ n n 2mu+a - 2 m\ 



n 



l+i/2 



i>0 



+ (-l)"/ 2 /i- n 2mv+a - 2 h/n l 



(20) 



i=0 



where the amplitude G m = G^ m 'B /m\ is related to the amplitude defined in equation fllffD . The 
scaling function prediction for the amplitudes G m is [3| 
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Figure 9. Plots of the estimates for some of the amplitude combinations G m B 



In figure |H1 we have plotted the resulting estimates for some of the amplitude ratios. Clearly the 
numerical results are fully consistent with the theoretical predictions. 



6. Summary and conclusions 

In this paper we have studied series for self-avoiding walks and polygons on the honeycomb 
lattice, including series for metric properties and moments of the area-distribution for polygons. 
We used various methods from Tony Guttmann's tool-kit to analyse the series. The connective 
constant, critical exponents and many amplitude combinations are known exactly, making it the 
perfect test-bed for series analysis techniques. 

In section 0] we used differential approximants to obtain estimates for the singularities and 
exponents of the SAP and SAW generating functions. Analysis of the SAP series (section 
14. lj) yielded very accurate estimates for the critical point x\ = 0.2928932186(5) and exponent 
2 — a = 1.5000004(10). The estimates agree with the conjectured exact values x 2 c = 1/(2 + \/2) = 
0.292893218813... and 2 — a = 3/2. In addition we found clear evidence of a non-physical 
singularity on the negative axis at x = X- = —0.412305(5) with an associated critical exponent 
a_ = 3/2. The analysis of the SAW series ( section I4.2JI also yielded estimates consistent with 
the predictions of the exact values. In this case there was a non-physical singularity at x = —x c 
as well as a pair of complex singularities at x = ±0.64215(15)i. So the excellent agreement 
is particularly impressive in light of the quite complicated functional form of the generating 
function which has at least three non-physical singularities as well as non-analytical corrections 
to scaling. So the walk series is obviously much harder to analyse and it is a great testament 
to the method of differential approximants that the analysis nevertheless yields such accurate 
estimates. 

In section we looked closer at the asymptotic form of the coefficients. In particular we 
obtained accurate estimates for the leading critical amplitudes A and B. Our method of analysis 
consisted in fitting the coefficients to an assumed asymptotic form. In section I5~T1 we analysed 
the SAP series and demonstrated clearly that in fitting to the coefficients we cannot ignore 
the singularity at x = X- even though it is exponentially suppressed asymptotically. After 
inclusion of this term estimates for the leading amplitude B were well behaved when including 
only analytic corrections and we found B = 1.2719299(1). We argued that this behaviour was 
consistent with a corrections-to-scaling exponent Ai being half-integer valued and in particular 
consistent with the prediction by Nienhuis that Ai = 3/2 (in section we showed the absence 
of a term with Ai = 11/16). In the analysis of the SAW series we discovered some subtleties 
regarding the inclusion of the terms arising from the complex pair of singularities. Despite 
a quite complicated asymptotic form (|16jl taking into account all the singularities and the 
corrections-to-scaling exponent Ai = 3/2 we could still obtain a quite accurate amplitude 
estimate A = 1.1449355(5). This analysis clearly shows that it is possible to probe quite deeply 



into the asymptotic behaviour of the series coefficients and in particular to distinguish between 
different corrections to scaling. 

E-mail or WWW retrieval of series 

The series for the generating functions studied in this paper can be obtained via e-mail by 
sending a request to I.Jensen@ms.unimelb.edu.au or via the world wide web on the URL 
http://www.ms.unimelb.edu.au/~iwan/ by following the instructions. 
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